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The recent formulation of the relativistic Thomas-Fermi model within the Feynman-Metropolis- 
Teller theory for compressed atoms is applied to the study of general relativistic white dwarf equilib- 
rium configurations. The equation of state, which takes into account the /3-equilibrium, the nuclear 
and the Coulomb interactions between the nuclei and the surrounding electrons, is obtained as a 
function of the compression by considering each atom constrained in a Wigner-Seitz cell. The con- 
tribution of quantum statistics, weak, nuclear, and electromagnetic interactions is obtained by the 
determination of the chemical potential of the Wigner-Seitz cell. The further contribution of the 
general relativistic equilibrium of white dwarf matter is expressed by the simple formula ^/goofJ-vs ~ 
constant, which links the chemical potential of the Wigner-Seitz cell /Xws with the general relativis- 
tic gravitational potential goo at each point of the configuration. The configuration outside each 
Wigner-Seitz cell is strictly neutral and therefore no global electric field is necessary to warranty 
the equilibrium of the white dwarf. These equations modify the ones used by Chandrasekhar by 
taking into due account the Coulomb interaction between the nuclei and the electrons as well as 
inverse /3-decay. They also generalize the work of Salpeter by considering a unified self-consistent 
approach to the Coulomb interaction in each Wigner-Seitz cell. The consequences on the numerical 
value of the Chandrasekhar-Landau mass limit as well as on the mass-radius relation of *He, ^^C, 
and ^®Fe white dwarfs are presented. All these effects should be taken into account in processes 
requiring a precision knowledge of the white dwarf parameters. 

Keywords: Relativistic Thomas-Fermi model - equation of state of white dwarf matter - white dwarf equi- 
librium configurations in general relativity 



I. INTRODUCTION 

The necessity of introducing the Fermi-Dirac statis- 
tics in order to overcome some conceptual difficulties 
in explaining the existence of white dwarfs leading to 
the concept of degenerate stars was first advanced by 
R. H. Fowler in a classic paper Following that work, 
E. C. Stoner 0] introduced the effect of special relativ- 
ity into the Fowler considerations and he discovered the 
critical mass of white dwarfs |56j] 
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where Mpi — ^Jhc/G ~ 10^^ g is the Planck mass, to„ 
is the neutron mass, and /i = A/Zsa2is the average 
molecular weight of matter which shows explicitly the 
dependence of the critical mass on the chemical compo- 
sition of the star. 

Following the Stoner's work, S. Chandrasekhar Q [s^l 
pointed out the relevance of describing white dwarfs by 
using an approach, initiated by E. A. Milne Q, of using 
the mathematical method of the solutions of the Lane- 
Emden polytropic equations [Hj. The same idea of using 
the Lane-Emden equations taking into account the spe- 
cial relativistic effects to the equilibrium of stellar matter 



for a degenerate system of fermions, came independently 
to L. D. Landau p. Both the Chandrasekhar and Lan- 
dau treatments were explicit in pointing out the existence 
of the critical mass 
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where the first numerical factor on the right hand 
side of Eq. ^ comes from the boundary condition 
~{r'^du/dr)r=R — 2.015 (see last entry of Table 7 on 
Pag. 80 in [5] ) of the n = 3 Lane-Emden polytropic equa- 
tion. Namely for M > M^j^"'^, no equilibrium configu- 
ration should exist. 

Landau rejected the idea of the existence of such a crit- 
ical mass as a "ridiculous tendency" [6]. Chandrasekhar 
was confronted by a lively dispute with A. Eddington on 
the basic theoretical assumptions he adopted [58] (see Q 
for details). 

Some of the basic assumptions adopted by Chan- 
drasekhar and Landau in their idealized approach were 
not justified e.g. the treatment of the electron as a free- 
gas without taking into due account the electromagnetic 
interactions, as well as the stability of the distribution of 
the nuclei against the gravitational interaction. It is not 
surprising that such an approach led to the criticisms of 
Eddington who had no confidence of the physical founda- 
tion of the Chandrasekhar work . It was unfortunate 
that the absence of interest of E. Fermi on the final evo- 
lution of stars did not allow Fermi himself to intervene in 
this contention and solve definitely these well-posed theo- 
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retical problems jsj. Indeed, we are showing in this article 
how the solution of the conceptual problems of the white 
dwarf models, left open for years, can be duly addressed 
by considering the relativistic Thomas-Fermi model of 
the compressed atom (see Subsec FlIEI and Sec. IIVI) . 

The original work on white dwarfs was motivated by 
astrophysics and found in astrophysics strong observa- 
tional support. The issue of the equilibrium of the elec- 
tron gas and the associated component of nuclei, taking 
into account the electromagnetic, the gravitational and 
the weak interactions is a theoretical physics problem, 
not yet formulated in a correct special and general rela- 
tivistic context. 

One of the earliest alternative approaches to 
the Chandrasekhar-Landau work was proposed by 
E. E. Salpeter in 1961 [9]. He followed an idea orig- 
inally proposed by Y. I. Frenkel [l3|: to adopt in the 
study of white dwarfs the concept of a Wigner-Seitz cell. 
Salpeter introduced to the lattice model of a point-like 
nucleus surrounded by a uniform cloud of electrons, cor- 
rections due to the non-uniformity of the electron distri- 
bution (see Subsec fll CI for details). In this way Salpeter 
obtained an analytic formula for the total energy in 
a Wigner-Seitz cell and derived the corresponding equa- 
tion of state of matter composed by such cells, pointing 
out explicitly the relevance of the Coulomb interaction. 

The consequences of the Coulomb interactions in the 
determination of the mass and radius of white dwarfs, 
was studied in a subsequent paper by T. Hamada and 
E. E. Salpeter [ll| by using the equation of state con- 
structed in P . They found that the critical mass of white 
dwarfs depends in a nontrivial way on the specific nuclear 
composition: the critical mass of Chandrasekhar-Landau 
which depends only on the mass to charge ratio of nuclei 
A/ Z , now depends also on the proton number Z . 

This fact can be seen from the approximate expression 
for the critical mass of white dwarfs obtained by Hamada 
and Salpeter [ll| in the ultrarelativistic limit for the elec- 
trons 
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being Ps the pressure of the Wigner-Seitz cell obtained 
by Salpeter in (see Subsec. IIIC|) and Pch is the pres- 
sure of a free-electron fluid used by Chandrasekhar (see 
Subsec. Ill A[) . The ratio Ps/Pch is a function of the 
number of protons Z (see Eq. (20) in Q) and it satis- 
fies Ps/Pch < 1- Consequently, the effective molecular 
weight satisfies /Zeff > /i and the critical mass of white 
dwarfs turns to be smaller than the original one obtained 
by Chandrasekhar-Landau (see Eq. ([2])). 

In the mean time, the problem of the equilibrium gas 
in a white dwarf taking into account possible global elec- 
tromagnetic interactions between the nucleus and the 



electrons was addressed by E. Olson and M. Bailyn in 
[l^ [Tsj . They well summarized the status of the prob- 
lem: ''^Traditional models for the white dwarf are non- 
relativistic and electrically neutral . . . although an electric 
field is needed to support the pressureless nuclei against 
gravitational collapse, the star is treated essentially in 
terms of only one charge component, where charge neu- 
trality is assumed " . Their solution to the problem in- 
vokes the breakdown of the local charge neutrality and 
the presence of an overall electric field as a consequence of 
treating also the nuclei inside the white dwarf as a fiuid. 
They treated the white dwarf matter through a two-fiuid 
model not enforcing local charge neutrality. The closure 
equation for the Einstein-Maxwell system of equations 
was there obtained from a minimization procedure of the 
mass-energy of the configuration. This work was the first 
pointing out the relevance of the Einstein-Maxwell equa- 
tions in the description of an astrophysical system by 
requiring global and non local charge neutrality. As we 
will show here, this interesting approach does not apply 
to the case of white dwarfs. It represents, however, a new 
development in the study of neutron stars (see e.g. [l^ ) 

An alternative approach to the Salpeter treatment of 
a compressed atom was reconsidered in |15| by applying 
for the first time to white dwarfs a relativistic Thomas- 
Fermi treatment of the compressed atom introducing a 
finite size nucleus within a phenomenological description 
(see also [Tg}). 

Recentl y, t he study of a compressed atom has been re- 
visited in [TtI by extending the global approach of Feyn- 
man. Metropolis and Teller 18| taking into account weak 
interactions. This treatment takes also into account all 
the Coulomb contributions duly expressed relativistically 
without the need of any piecewise description. The rela- 
tivistic Thomas-Fermi model has been solved by impos- 
ing in addition to the electromagnetic interaction also the 
weak equilibrium between neutrons, protons and elec- 
trons self-consistently. This presents some conceptual 
differences with respect to previous approaches and can 
be used in order both to validate and to establish their 
limitations. 

In this article we apply the considerations presented in 
[l7| of a compressed atom in a Wigner-Seitz cell to the 
description of non-rotating white dwarfs in general rel- 
ativity. This approach improves all previous treatments 
in the following aspects: 

1. In order to warranty self-consistency with a rela- 
tivistic treatment of the electrons, the point-like as- 
sumption of the nucleus is abandoned introducing a 
finite sized nucleus |l7[ • We assume for the mass as 



well as for charge to mass ratio of the nucleus their 
experimental values instead of using phenomeno- 
logical descriptions based on the semi-empirical 
mass-formula of Weizsacker (see e.g. (isl [l6|). 

2. The electron-electron and electron- nucleus Cou- 
lomb interaction energy is calculated without any 
approximation by solving numerically the relativis- 
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tic Thomas-Fermi equation for selected energy- 
densities of the system and for each given nuclear 
composition. 

3. The energy-density of the system is calculated tak- 
ing into account the contributions of the nuclei, of 
the Coulomb interactions as well as of the relativis- 
tic electrons; the latter being neglected in all previ- 
ous treatments. This particular contribution turns 
to be very important at high-densities and in par- 
ticular for light nuclear compositions e.g. ^Hc and 

4. The /3-equilibrium between neutrons, protons, and 
electrons is also taken into account leading to a self- 
consistent calculation of the threshold density for 
triggering the inverse /3-decay of a given nucleus. 

5. The structure of the white dwarf configurations is 
obtained by integrating the general relativity equa- 
tions of equilibrium. 

6. Due to 4) and 5) we are able to determine if the 
instability point leading to a maximum stable mass 
of the non-rotating white dwarf is induced by the 
inverse /3-decay instability of the composing nuclei 
or by general relativistic effects. 

Paradoxically, after all this procedure which takes into 
account many additional theoretical features generaliz- 
ing the Chandrasekhar-Landau and the Hamada and 
Salpeter works, a most simple equation is found to be 
fulfilled by the equilibrium configuration in a spherically 
symmetric metric. Assuming the metric 



e-i-)c^dt' - e^^-Ur'' - r^dd^ ~ sin^ ddip^ , (5) 



we demonstrate how the entire system of equations de- 
scribing the equilibrium of white dwarfs, taking into ac- 
count the weak, the electromagnetic and the gravita- 
tional interactions as well as quantum statistics all ex- 
pressed consistently in a general relativistic approach, is 
simply given by 



/goo^ws = e''^''^^^fiws{r) = constant. 
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which links the chemical potential of the Wigner-Seitz 
cell ^ws: duly solved by considering the relativistic 
Feynman- Metropolis- Teller model following to the 
general relativistic gravitational potential at each point 
of the configuration. The overall system outside each 
Wigner-Seitz cell is strictly neutral and no global electric 
field exists, contrary to the results reported in [1^. The 
same procedure will apply as well to the case of neutron 
star crusts. 

The article is organized as follows. In Sec.|ll]we sum- 
marize the most common approaches used for the descrip- 
tion of white dwarfs and neutron star crusts: the uniform 
approximation for the electron fiuid (see e.g. @); the 
often called lattice model assuming a point-like nucleus 
surrounded by a uniform electron cloud (see e.g. [l9|): 



the generalization of the lattice model due to Salpeter 
@; the Feynman, Metropolis and Teller approach 
based on the the non-relativistic Thomas-Fermi model 
of compressed atoms and, the relativistic generalization 
of the Feynman- Metropolis- Teller treatment recently for- 
mulated in [TtI . 

In Sec. mil we formulate the general relativistic equa- 
tions of equilibrium of the system and show how, from 
the self-consistent definition of chemical potential of the 
Wigner-Seitz cell and the Einstein equations, comes the 
equilibrium condition given by Eq. ([6]). In addition, we 
obtain the Newtonian and the first-order post-Newtonian 
equations of equilibrium. 

Finally, we show in Sec. IIVI the new results of the nu- 
merical integration of the general relativistic equations of 
equilibrium and discuss the corrections to the Stoner crit- 
ical mass M^j.*°"°'', to the Chandrasekhar-Landau mass 
limit M^f}.~^, as well as to the one of Hamada and 
Salpeter M^f^^, obtained when all interactions are fully 
taken into account through the relativistic Feynman- 
Metropolis- Teller equation of state [l7| . 



II. THE EQUATION OF STATE 

There exists a large variety of approaches to model 
the equation of state of white dwarf matter, each one 
characterized by a different way of treating or neglecting 
the Coulomb interaction inside each Wigner-Seitz cell, 
which we will briefiy review here. Particular attention 
is given to the calculation of the self-consistent chemical 
potential of the Wigner-Seitz cell /iws, which plays a very 
important role in the conservation law ^ that we will 
derive in Sec. IIIII 



A. The uniform approximation 

In the uniform approximation used by Chandrasekhar 
3, the electron distribution as well as the nucleons are 
assumed to be locally constant and therefore the condi- 
tion of local charge neutrality 
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where Ar is the average atomic weight of the nucleus, 
is applied. Here n^r denotes the nucleon number den- 
sity and Z is the number of protons of the nucleus. 
The electrons are considered as a fully degenerate free- 
gas and then described by Fermi-Dirac statistics. Thus, 
their number density Ue is related to the electron Fermi- 
momentum Pf by 
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and the total electron energy-density and electron pres- 
sure are given by 
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where we have introduced the dimensionless Fermi mo- 
mentum Xe — Pf /(wec) with iTie the electron rest- mass. 

The kinetic energy of nucleons is neglected and there- 
fore the pressure is assumed to be only due to electrons. 
Thus the equation of state can be written as 

■^^unif = £n +£e~ -^M„C^ne -|- £e , (H) 
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where M,j = 1.6604 x 10~24 g jg the unified atomic mass 
and £e and Pg are given by Eqs. ([9t- p0|) . 

Within this approximation, the total self-consistent 
chemical potential is given by 
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is the electron free-chemical potential. 

As a consequence of this effective approach which does 
not take into any account the Coulomb interaction, it 
is obtained an effective one-component electron-nucleon 
fluid approach where the kinetic pressure is given by 
electrons of mass nie and their gravitational contribu- 
tion is given by an effective mass {Ar/Z)Mu attached to 
each electron (see e.g. [13] )• This is even more evident 
when the electron contribution to the energy-density in 
Eq. ([Tl]) is neglected and therefore the energy-density is 
attributed only to the nuclei. Within this approach fol- 
lowed by Chandrasekhar , the equation of state reduces 
to 
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B. The lattice model 



The first correction to the above uniform model, cor- 
responds to abandon the assumption of the electron- 
nucleon fluid through the so-called "lattice" model which 



introduces the concept of Wigner-Seitz cell: each cell con- 
tains a point-like nucleus of charge +Ze with A nucleons 
surrounded by a uniformly distributed cloud of Z fully- 
degenerate electrons. The global neutrality of the cell is 
guaranteed by the condition 



Z = Kvsrie 



(17) 



where riws = 1/Krs is the Wigner-Seitz cell density and 
Vws = 47ri?^g/3 is the cell volume. 

The total energy of the Wigner-Seitz cell is modified 
by the inclusion of the Coulomb energy, i.e 
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where £unif is given by Eq. (|lll) and E^-n and Ee-e are 
the electron-nucleus and the electron-electron Coulomb 
energies 
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5 Pws 

The self-consistent pressure of the Wigner-Seitz cell is 
then given by 

P _ lEc 

Pl - -7^7^ - Punif + , [2Z) 

where Punif is given by Eq. ([T^. It is worth to recall 
that the point-like assumption of the nucleus is incom- 
patible with a relativistic treatment of the degenerate 
electron fluid (see [HI [13] for details). Such an inconsis- 
tency has been traditionally ignored by applying, within 
a point-like nucleus model, the relativistic formulas Q 
and (|10p and their corresponding ultrarelativistic limits 
(see e.g. 9]). 

The Wigner-Seitz cell chemical potential is in this case 
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By comparing Eqs. (|T2|) and (|22l) we can see that the 
inclusion of the Coulomb interaction results in a decreas- 
ing of the pressure of the cell due to the negative lat- 
tice energy Ec- The same conclusion is achieved for the 
chemical potential from Eqs. (fT3|) and (|23l) . 



C. Salpeter approach 

A further development to the lattice model came from 
Salpeter Q whom studied the corrections due to the non- 
uniformity of the electron distribution inside a Wigner- 
Seitz cell. 
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Following the Chandrasekhar |3| approximation, 
Salpeter also neglects the electron contribution to the 
energy-density. Thus, the first term in the Salpeter for- 
mula for the energy of the cell comes from the nuclei 
energy (|15p . The second contribution is given by the 
Coulomb energy of the lattice model (fTQ]). The third 
contribution is obtained as follows: the electron density 
is assumed as ne[l -|- e(r)], where rte = 3Z/(47ri?^g) is 
the average electron density as given by Eq. p7)) . and 
e(r) is considered infinitesimal. The Coulomb potential 
energy is assumed to be the one of the point-like nucleus 
surrounded by a uniform distribution of electrons, so the 
correction given by e(r) on the Coulomb potential is ne- 
glected. The electron distribution is then calculated at 
first-order by expanding the relativistic electron kinetic 
energy 



- mec^, (24) 
about its value in the uniform approximation 



= ^;i2c2(3^2„^)2/3+^2c4 _ ^^^2 ^ (25) 

considering as infinitesimal the ratio eV/E^ between the 
Coulomb potential energy eV and the electron Fermi en- 
ergy 



= y [cPf (r)]2 + m2c4 - m^c^ - eV . (26) 

The influence of the Dirac electron-exchange correc- 
tion [23'] on the equation of state was also considered by 
Salpeter [9,]. However, adopting the general approach of 
Migdal et al. [2^, it has been shown that these effects 
are negligible in the relativistic regime flT'l . We will then 
consider here only the major correction of the Salpeter 
treatment. 

The total energy of the Wigner-Seitz cell is then given 
by (see for details) 
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where Ech = ^chK/s, Ec is given by Eq. (HH), is 
given by Eq. and a = e"^ / (he) is the fine structure 

constant. 

Correspondingly, the self-consistent pressure of the 
Wigner-Seitz cell is 



Ps = Pl+ pSf , 



where 
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The Wigner-Seitz cell chemical potential can be then 
written as 
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From Eqs. ((29)) and ((3T|) . we see that the inclusion of 
each additional Coulomb correction results in a further 
decreasing of the pressure and of the chemical poten- 
tial of the cell. The Salpeter approach is very interest- 
ing in identifying piecewise Coulomb contribution to the 
total energy, to the total pressure and, to the Wigner- 
Seitz chemical potential. However, it does not have the 
full consistency of the global solutions obtained with the 
Feynman-Metropolis- Teller approach Is'] and its gener- 
alization to relativistic regimes ^17i] which we will discuss 
in detail below. 



D. The Feynman-Metropolis- Teller treatment 

Feynman, Metropolis and Teller [l^ showed how to 
derive the equation of state of matter at high pressures by 
considering a Thomas-Fermi model confined in a Wigner- 
Seitz cell of radius i?ws- 

The Thomas-Fermi equilibrium condition for degener- 
ate non-relativistic electrons in the cell is expressed by 



{Pi? 
2me 



E, 
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(32) 



where V denotes the Coulomb potential and E^ denotes 
the Fermi energy of electrons, which is positive for con- 
figurations subjected to external pressure, namely, for 
compressed cells. 

Defining the function 0(r) by eV{r)+E^ = e'^ Z(f>{r)/r, 
and introducing the dimensionless radial coordinate 77 
by r = b-q, where b = {37r)^^^{\e/a)2~''^^Z-^/^, being 
Ae = h/{mec) the electron Compton wavelength; the 
Poisson equation from which the Coulomb potential V 
is calculated self-consistently becomes 



(33) 



The boundary conditions for Eq. (j33l) follow from the 
point-like structure of the nucleus 0(0) = 1 and, from 
the global neutrality of the Wigner-Seitz cell </'(??o) = 
riodcj) / dri\ri=r)a , where ryo defines the dimensionless radius 
of the Wigner-Seitz cell by 770 — R^is/b. 

For each value of the compression, e.g. 770, it corre- 
sponds a value of the electron Fermi energy E^ and 
a different solution of Eq. which determines the 

self-consistent Coulomb potential energy eV as well as 
the self-consistent electron distribution inside the cell 
through 



3/2 



(34) 
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In the non-relativistic Thomas-Fermi model, the total 
energy of the Wigner-Seitz cell is given by (see [H, [11] 
for details) 
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where A1n{Z,A) is the nucleus mass, Sein-eir)] is given 
by Eq. © and E^-n and E^-e are the electron-nucleus 
Coulomb energy and the electron-electron Coulomb en- 
ergy, which are given by 
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From Eqs. ([57]) and we recover the well-known 
relation between the total kinetic energy and the total 
Coulomb energy in the Thomas- Fermi model [3 HEl 

4^)=i?r'[ne(i?w.)]-^i?c, (41) 

where E'^"'^ [ne(i?„s)] is the non-relativistic kinetic energy 
of a uniform electron distribution of density ne(-Rws), i-e. 



i?r'K(i?ws)] = -^>e(i?ws) 
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with Z* defined by 



Z* — Vft,s"-e(-Rws) 7 



(42) 



(43) 



and Me(i?ws) = ?l2[37r2ne(i?ws)]^/V(2TOe). 

The self-consistent pressure of the Wigner-Scitz cell 
give n by the non-relativistic Thomas-Fermi model is (see 
[il mi for details) 
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The pressure of the Thomas-Fermi model (|44)) is 
equal to the pressure of a free-electron distribution 
of density ne(i?ws)- Being the electron density in- 
side the cell a decreasing function of the distance from 



the nucleus, the electron density at the cell boundary, 
^e(-Rws)j is smaller than the average electron distribu- 
tion 3Z/(47ri?^g). Then, the pressure given by (j44|) is 
smaller than the one given by the non-relativistic version 
of Eq. ([T0|) of the uniform model of Subsec. Ill Al Such a 
smaller pressure, although faintfully given by the expres- 
sion of a free-electron gas, contains in a self-consistent 
fashion all the Coulomb effects inside the Wigner-Seitz 
cell. 

The chemical potential of the Wigner-Seitz cell of the 
non-relativistic Thomas- Fermi model can be then written 
as 

MTF = Mn{Z, A)^ + Z*ii,{R^s) + \ec , (45) 

where we have used Eqs. (PT|) - (|15)) . 

Integrating by parts the total number of electrons 



Z = 



4Trr^n,{r)dr = Z* + I{R^,) , (46) 



where 
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3 dr ' 



(47) 



we can rewrite finally the following semi-analytical ex- 
pression of the chemical potential (|45|) of the cell 
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where /x""'^ is the electron free-chemical potential (fT4|) 
calculated with the average electron density, namely, the 
electron chemical potential of the uniform approxima- 
tion. The function /(i?ws) depends explicitly on the gra- 
dient of the electron density, i.e. on the non-uniformity 
of the electron distribution. 

In the limit of absence of Coulomb interaction both 
the last term and the function /(i?„s) in Eq. (|48l) vanish 
and therefore in this limit /itf reduces to 



MTF — >■ Munif , 



(49) 



where /^unif is the chemical potential in the uniform ap- 
proximation given by Eq. ([Tc 



E. The relativistic Feynman-Metropolis-Teller 
treatment 

We recall now how the above classic Feynman, 
Metropolis, and Teller treatment of compressed atoms 
has been recently generalized to relativistic regimes (see 
(Tt] for details). One of the main differences in the rel- 
ativistic generalization of the Thomas-Fermi equation is 
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that, the point-like approximation of the nucleus, must 
be abandoned since the relativistic equilibrium condition 
of compressed atoms 



constant > , 
(50) 

would lead to a non-integrable expression for the electron 
density near the origin (see e.g. [2T!, 125]). 

It is then assumed a constant distribution of protons 
confined in a radius Rc defined by 



1/3 



(51) 



where Att ~ h/{m.„c) is the pion Compton wavelength. 
If the system is at nuclear density A w {r()/\-,;){A/ZY/'^ 
with ro ~ 1.2 fm. Thus, in the case of ordinary nuclei 
(i.e., for A/Z « 2) we have A « 1. Consequently, the 
proton density can be written as 



np{r) 



1 



(r^Rc), (52) 



where 9{r — Rc) denotes the Heaviside function centered 
at Rc- The electron density can be written as 



neir) 



F\3 



1 



(Pi) 



3/2 



(53) 

where V — eV + and we have used Eq. (|50|) . 

The overall Coulomb potential satisfies the Poisson 
equation 



V'^V{r) = -47re [np[r) - n^ir)] , 



(54) 



with the boundary conditions dV / dr\r=R^^ = and 
V{R^s) = due to the global charge neutrality of the 
cell. 

By introducing the dimensionless quantities x — r /X^^^ 

Xc = Rc/XjT, x/f = ^('')/(^c) and replacing the particle 
densities (|52|) and ((53)) into the Poisson equation (|54|. it 
is obtained the relativistic Thomas-Fermi equation [26| 



1 d^xjx) 
3x dx"^ 



4Q! 

97r 



X^{x) me x{x) 



3/2 



(55) 



which must be integrated subjected to the boundary 
conditions x(0) = 0, xi^vis) > and dx/ dx\x=x^s — 
X{x^s)/xws, where a;ws = i?ws/A^. 

The neutron density n„(r), related to the neutron 
Fermi momentum P.^ — (37r2/i^n„)^/'^, is determined by 
imposing the condition of beta equilibrium 



mpC^ + eV{r) + E^ 



(56) 



Relativistic FMT treatment 
Uniform approximation 




250 



300 



FIG. 1: The electron number density Tie in units of the average 
electron number density no = 3Z/{4nRl,s) inside a Wigner- 
Seitz cell of ^'^C. The dimensionless radial coordinate is a; = 
r/A,r and Wigner-Seitz cell radius ~ 255 corresponding 

to a density of ~ 10* g/cm^. The solid curve corresponds 
to the relativistic Feynman-Metropolis- Teller treatment and 
the dashed curve to the uniform approximation. The electron 
distribution for different levels of compression as well as for 
different nuclear compositions can be found in [l7j |. 



subjected to the baryon number conservation equation 



A 



47rr^[np(r) + n„(r)]c?r . 



(57) 



In Fig. [T] we see how the relativistic generalization of 
the Feynman-Metropolis- Teller treatment leads to elec- 
tron density distributions markedly different from the 
constant electron density approximation. The electron 
distribution is far from being uniform as a result of the 
solution of Eq. ([55]) , which takes into account the electro- 
magnetic interaction between electrons and between the 
electrons and the finite sized nucleus. Additional details 
are given in [13]. 

V. S. Popov et al. 27[ have shown how the solu- 
tion of the relativistic Thomas-Fermi equation (|55|) to- 



gether with the self-consistent implementation of the /3- 
equilibrium condition (|56l) leads, in the case of zero elec- 
tron Fermi energy {E^ = 0), to a theoretical prediction of 
the /3-equilibrium line, namely a theoretical Z-A relation. 
Within this model the mass to charge ratio A/Z of nuclei 
is overestimated, e.g. in the case of *He the overestimate 
is ~ 3.8%, for ~ 7.9%, for ^^O - 9.52%, and for ^^Fe 
~ 13.2%. These discrepancies are corrected when the 
model of the nucleus considered above is improved by ex- 
plicitly including the effects of strong interactions. This 
model, however, illustrates how a self-consistent calcula- 
tion of compressed nuclear matter can be done including 
electromagnetic, weak, strong as well as special relativis- 
tic effects without any approximation. This approach 
promises to be useful when theoretical predictions are 
essential, for example in the description of nuclear mat- 
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ter at very high densities, e.g. nuclei close and beyond 
the neutron drip line. 

The densities in white dwarf interiors are not highly 
enough to require such theoretical predictions. There- 
fore, in order to ensure the accuracy of our results we use 
for {Z, A) , needed to solve the relativistic Thomas- Fermi 
equation (|55|) . as well as for the nucleus mass Mm{Z, A), 
their known experimental values. In this way we take 
into account all the effects of the nuclear interaction. 

Thus, the total energy of the Wigner-Seitz cell in the 
present case can be written as 



-C'FMT ~ 



Er 



being 



En 



E. 



(e) 



MNiZ,A)c^, 

47rr^ {8e — mf,nf,)dr , 



(58) 

(59) 
(60) 



Ec 



1 



4:TTr'^ e[np{r) ~ ne{r)]V{r)dr , (61) 



where Mn{Z,A) = ArMu is the experimental nucleus 
mass, e.g. for ''He, ^^C, ^^O and ^^Fe we have Ar = 
4.003, 12.01, 16.00 and 55.84 respectively. In Eq. §^ 
the integral is evaluated only outside the nucleus (i.e. 
for r > i?c) in order to avoid a double counting with the 
Coulomb energy of the nucleus already taken into account 
in the nucleus mass ((59|) . In order to avoid another dou- 
ble counting we subtract to the electron energy-density 
£e in Eq. (j60p the rest-energy density nieC^ne which is 
also taken into account in the nucleus mass (j59|) . 

The total pressure of the Wigner-Seitz cell is given by 



FMT 



(62) 



where Pe[ne{Rws)] is the relativistic pressure pO)) com- 
puted with the value of the electron density at the bound- 
ary of the cell. 

The electron density at the boundary R„s in the rel- 
ativistic Feynman-Metropolis- Teller treatment is smaller 
with respect to the one given by the uniform density ap- 
proximation (see Fig. [T|). Thus, the relativistic pressure 
(l62|l gives systematically smaller values with respect to 
the uniform approximation pressure (jlOp as well as with 
respect to the Salpeter pressure (j^ . 

In Fig. [5] we show the ratio between the relativistic 
Feynman-Metropolis- Teller pressure -Pp^T f62)) and the 
Chandrasekhar pressure Pch pH]) and the Salpeter pres- 
sure Ps (HH) in the case of ^^C. It can be seen how -Pp^T 
is smaller than Pch for all densities as a consequence of 
the Coulomb interaction. With respect to the Salpeter 
case, we have that the ratio Pfmt/^s approaches unity 
from below at large densities as one should expect. 

However, at low densities < IC'-IO^ g/cm'^, the ra- 
tio becomes larger than unity due to the defect of 
the Salpeter treatment which, in the low density non- 
relativistic regime, leads to a drastic decrease of the pres- 
sure and even to negative pressures at densities < 10^ 



1.05 



1.00 



0.95 



0.90 



-fFMT/-fcii 



10-' 



10" 



10-' 



10" 

P (g/c 



10' 



10" 



10' 



FIG. 2: Ratio of the pressures in the different treatments as a 
function of the density for ^^C white dwarfs (see Table|lJ. The 
solid curve corresponds to the ratio between the relativistic 
Feynman-Metropolis- Teller pressure Pfmt given by Eq. (|62p 



and the Chandrasekhar pressure Pch given by Eq. (|10|) . The 
dashed curve corresponds to the ratio between the relativistic 
Feynman-Metropolis- Teller pressure Pfmt given by Eq. (|62p 
and the Salpeter pressure Ps given by Eq. (|29p . 



Ch 



Ps 



prcl 
-'FMT 



10 
40 
70 
10^ 
10^ 

lo-* 

10^ 
10^ 
10^ 
10« 
10« 
10^° 



1.46731 X 
1.47872 X 
3.75748 X 
6.80802 X 
3.15435 X 
1.45213 X 
6.50010 X 
2.62761 X 
8.46101 X 
2.15111 X 
4.86236 X 
1.05977 X 



lO^-* 
lO^s 
10^5 
10^5 
10^^ 
lO^^ 
10^0 
10^ 
10^ 



|22 
10^6 



-1.35282 X 10^^ 
4.60243 X 10" 
1.60860 X 10^^ 
3.34940 X 10^^ 
2.40646 X 10^^ 
1.28976 X 10^^ 
6.14494 X 10^° 
2.54932 X 10^^ 
8.28899 X 10^^ 
2.11375 X 10^^ 
4.78170 X 10^^ 
1.04239 X 10^** 



4.54920 X 10" 
7.09818 X 10" 
2.05197 X 10^^ 
3.90006 X 10^^ 
2.44206 X lO^'' 
1.28965 X 10^" 
6.13369 X 10^° 
2.54431 X 10^^ 
8.27285 X 10^^ 
2.10896 X 10^^^ 
4.76613 X 10^® 
1.03668 X 10^* 



TABLE I: Equation of state for ^^C within the different 
treatments. The pressure in the uniform approximation for 
fi — 2 is Pch, the Salpeter pressure is Ps and the relativistic 
Feynman-Metropolis- Teller pressure is Pfmt- The units for 
the density are g/cm^ and for the pressure dyn/cm^. 



g/cm^ or higher for heavier nuclear compositions e.g. ^^Fe 
(see [i, 113 and Table H]). This is in contrast with the 
relativistic Feynman-Metropolis- Teller treatment which 
matches smoothly the classic Feynman-Metropolis- Teller 
equation of state in that regime (see [TtI] for details). 

No analytic expression of the Wigner-Seitz cell chem- 
ical potential can be given in this case, so we only write 
its general expression 



,,rol prcl 

PFMT — -'^FMT ■ 



prcl Y 



(63) 
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where -Ep^x ^'^'^ -^fmt ^^'^ given by Eqs. ([58]) and ((62)) 
respectively. The above equation, contrary to the non- 
relativistic formula (|45p . in no way can be simplified in 
terms of its uniform counterparts. However, it is easy 
to check that, in the limit of no Coulomb interaction 
ne{Rws) 3Z/(47ri?3j, Ec 0, and Ek fchK-s 
and, neglecting the nuclear binding and the proton- 
neutron mass difference, we finally obtain 



Mfmt 



(64) 



as it should be expected. 

Now we summarize how the equation of state of com- 
pressed nuclear matter can be computed in the Salpeter 
case and in the relativistic Feynman-Metropolis- Teller 
case, parameterized by the total density of the system: 

(i) For a given radius i?ws of the Wigner-Seitz cell the 
relativistic Thomas-Fermi equation (|55l) is integrated nu- 
merically and the density of the configuration is com- 
puted as p = Epl^rj,/[c'^V„s) where iJpMT the energy 



of the cell given by Eq. ([S? 

(ii) For that value of the density, the radius of the 
Wigner-Seitz cell in the Salpeter treatment is 



1/3 



(65) 



where Eq. ()15p has been used. On the contrary, in the 
relativistic Feynman-Metropolis- Teller treatment no an- 
alytic expression relating Wigner-Seitz cell radius and 
density can be written. 

(iii) From this Wigner-Seitz cell radius, or equivalently 
using the value of the density, the electron density in 
the Salpeter model is computed from the assumption of 
uniform electron distribution and the charge neutrality 
condition, i.e. Eq. (I15|) . In the relativistic Feynman- 
Metropolis- Teller treatment, the electron number den- 
sity at the boundary of the Wigner-Seitz cell is, following 
Eq. given by 



olFMT 



1 



37r2A3 



ws J 



+ 2 



3/2 



(66) 



where the function x{^) is the solution of the relativistic 
Thomas- Fermi equation (|55p . 

(iv) Finally, with the knowledge of the electron density 
at i?ws, the pressure can be calculated. In the Salpeter 
approach it is given by Eq. while in the relativistic 
Feynman-Metropolis- Teller case it is given by Eq. (p^ . 



III. GENERAL RELATIVISTIC EQUATIONS OF 
EQUILIBRIUM 

Outside each Wigner-Seitz cell the system is electri- 
cally neutral, thus no overall electric field exists. There- 
fore, the above equation of state can be used to calcu- 
late the structure of the star through the Einstein equa- 
tions. Introducing the spherically symmetric metric ([5]), 



the Einstein equations can be written in the Tolman- 
Oppenheimer-Volkoff form [H H 



dv(r) 
dr 

dM{r) 

dr 
dP{r) 
dr 



2G 4Trr^P{r)/c^ + M{r) 



47rr2— V- 
1 di>(r^ 



2GM(r) 



2 dr 



■[f(r)+P(r)] 



(67) 

(68) 
(69) 



where we have introduced the mass enclosed at the dis- 
tance r through e^^''^ = 1 - 2GM[r)/ {c^r), £{r) is the 
energy-density and P{r) is the total pressure. 

We turn now to demonstrate how, from Eq. (|69p . it 
follows the general relativistic equation of equilibrium 
dU, for the self-consistent Wigner-Seitz chemical poten- 
tial /iws- The first law of thermodynamics for a zero 
temperature fluid of N particles, total energy i?, total 
volume V , total pressure P = —dE/dV, and chemical 
potential /i = dE/dN reads 



dE = -PdV + ndN , 



(70) 



where the differentials denote arbitrary but simultaneous 
changes in the variables. Since for a system whose surface 
energy can be neglected with respect to volume energy, 
the total energy per particle E/N depends only on the 
particle density n = N/V, we can assume E/N as an ho- 
mogeneous function of first-order in the variables N and 
V and hence, it follows the well-known thermodynamic 
relation 

E = -PV + ^iN. (71) 
In the case of the Wigner-Seitz cells, Eq. ([TTt reads 



-P V 



(72) 



where we have introduced the fact that the Wigner-Seitz 
cells are the building blocks of the configuration and 
therefore we must put in Eq. (j7ip N^^s — 1. Through 
the entire article we have used Eq. (17^ to obtain from 
the knowns energy and pressure, the Wigner-Seitz cell 
chemical potential (see e.g. Eqs. ([T^ and (|23I) V From 
Eqs. (j70p and dTlT) we obtain the so-called Gibbs-Duhem 
relation 



dP = nd^ . 



(73) 



In a white dwarf the pressure P and the chemical po- 
tential fi are decreasing functions of the distance from 
the origin. Thus, the differentials in the above equations 
can be assumed as the gradients of the variables which, 
in the present spherically symmetric case, become just 
derivatives with respect to the radial coordinate r. From 
Eq. ([75)1 it follows the relation 



dP^ 
dr 



dr 



(74) 
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From Eqs. dSH), and we obtain 

n-wsir) ^^^!"^^^ = T~^«ws(r)/iws(r) , (75) 



dr 



2 dr 



which can be straightforwardly integrated to obtain the 
first integral 



g'^('')/2^^_^^^-) — constant. 



(76) 



The above equilibrium condition is general and it also 
applies^ for non-zero temperature configurations ( see 
e.g. 30]). In such a case, it can be shown that in ad- 
dition to the equilibrium condition (j76p the temperature 
of the system satisfies the Tolman isothermality condi- 
tion e''('')/2r(r) = constant [HI, HI]. 



A. The weak-field non-relativistic limit 

In the weak-field limit we have e'^^'^ « 1 -I- <&, where 
the Newtonian gravitational potential has been defined 
by $(r) = v{r)/2. In the non-relativistic mechanics limit 
c — >■ cxD, the chemical potential /iws — >■ /iws+AfwsC^, where 
/2ws denotes the non-relativistic free-chemical potential 
of the Wigner-Seitz cell and M„s is the rest-mass of the 
Wigner-Seitz cell, namely, the rest-mass of the nucleus 
plus the rest-mass of the electrons. Applying these con- 
siderations to Eq. (fTS)) we obtain 



Afws$ — constant . 



(77) 



Absorbing the Wigner-Seitz rest-mass energy A/^sC^ in 
the constant on the right-hand-side we obtain 



Aiws'I' = constant . 



(78) 



In the weak- field non-relativistic limit, the Einstein 
equations ([67)) - l|69|) reduce to 



d$(r) 

dr 
dM{r) 

dr 
dP{r) 

dr 



GM[r) 

47rr^p(r) 
GM{r) 



(79) 
(80) 
(81) 



where p{r) denotes the rest-mass density. The Eqs. (|79p - 
(150)) can be combined to obtain the gravitational Poisson 
equation 



d2$(r) 2 d<^{r) 



dr"^ 



dr 



AirGpir) . 



(82) 



In the uniform approximation (see Subsec. Ill A[) . the 
equilibrium condition given by Eq. ([75)) reads 



-M„$ = constant . 



(83) 



equation by the total number of electrons Z. This equi- 
librium equation is the classical condition of thermo- 
dynamic equilibrium assumed for non-relativistic white 
dwarf models (see e.g. \2d\ for details). 

Introducing the above equilibrium condition (|83p into 
Eq. (15^ . and using the relation between the non- 
relativistic electron chemical potential and the particle 
density = (2me)^/^/ie^^/(37r^fi,'^), we obtain 

d^fleir) 2dile{r) _ 'I'l'^nC'J'^^ArlZfra^G^^ 



dr"^ 



dr 



37rft3 



... 

which is the correct equation governing the equilibrium of 
white dwarfs within Newtonian gravitational theory pol | . 
It is remarkable that the equation of equilibrium ()84l) . 
obtained from the correct application of the Newtonian 
limit, does not coincide with the equation given by [3l.l33l- 
|35| , which, as correctly pointed out by , is a mixture of 
both relativistic and non-relativistic approaches. Indeed, 
the consistent relativistic equations should be Eq. ([76)) . 
Therefore a dual relativistic and non-relativistic equation 
of state was used by Chandrasekhar. The pressure on 
the left-hand-side of Eq. ([STt is taken to be given by 
relativistic electrons while, the term on the right-hand- 
side of Eq. ([80]) and ^ (or the source of Eq. ([82])), is 
taken to be the rest-mass density of the system instead of 
the total relativistic energy-density. Such a procedure is 
equivalent to take the chemical potential in Eq. ([78| as a 
relativistic quantity. As we have seen, this is inconsistent 
with the weak-field non-relativistic limit of the general 
relativistic equations. 



B. The Post-Newtonian limit 

Although quantitatively justifiable (see next section), 
the Chandrasekhar approach was strongly criticized by 
Eddington because it was conceptually unjustified. In- 
deed, if one were to treat the problem of white dwarfs 
approximately without going to the sophistications of 
general relativity, but including the effects of relativistic 
mechanics, one should use at least the equations in the 
post-Newtonian limit. The first-order post-Newtonian 
expansion of the Einstein equations ()67l) - ()69p in powers 
of P/£ and GM/{c^r) leads to the equilibrium equations 

M 



d4>(r) 

dr 
dM(r) 

dr 
dP{r) 

dr 



1 



£{r) 



1 



P{r) 



£{r) 



■Ar) 



GM{r) £{r) 
2GM{r)' 



dP{r) 
dr 



P{r) Airr^Pir) 
Jirj ^ M(r)c2 



(85) 
(86) 

) 

(87) 



where we have neglected the electron rest-mass with re- 
spect to the nucleus rest-mass and we have divided the 



where Eq. ([57)) is the post-Newtonian version of the 
Tolman-Oppenheimer-Volkoff equation ([M)) . 
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Replacing Eq. (|74)) into Eq. ([85]) we obtain 
P(r) 



1 - 



Sir) 



dfi^r) ^ £{r)/c^d^r) 
dr 



dr n^s{r) 

It is convenient to split the energy-density as £ = c? p+U 
where p = Af^s^-ws is the rest-energy density and U the 
internal energy-density. Thus, Eq 





Decay 


4 


^(.rclFMT 
Pcrit 


/3,unif 
Pcrit 


''He - 


► ^ H + n 4n 


20.596 


1.39 X 10" 


1.37 X 10" 


12^ 


-^i^B ^i2Be 


13.370 


3.97 X 10^° 


3.88 X 10^" 


16q 




10.419 


1.94 X 10^° 


1.89 X 10^" 






3.695 


1.18 X 10^ 


1.14 X 10^ 



dr 



becomes 



d$(r) P{r)dp^,{r) 



dr 



U/c^ d$(r) 
nws(r) dr 



£{r) 
= 0, 



dr 



TABLE II: Onset of inverse beta decay instability for ''He, 
^^C, ^''O and ^"^Fe. The experimental inverse /3-decay ener- 
gies are given in MeV and they have been taken from Table 
1 of [431. The corresponding critical density for the uniform 

given by Eq. (|90p . is given in 



electron density model, pfrf""'' 



(89) g/cm^ as well as the critical density pj? 



which is the differential post-Newtonian version of the 
equilibrium equation (|76l) and where the post-Newtonian 
corrections of equilibrium can be clearly seen. Applying 
the non-relativistic limit c — ?> oo to Eq. (|89l) : P/5 — > 



for the relativis- 
tic Feynman-Metropolis- Teller case. The numerical values of 
are taken from [31) see also [i^] 



0, U/(T -> 0, and /iws Mws' 
Newtonian equation of equilibrium ([75]) . 



IV. MASS AND RADIUS OF GENERAL 
RELATIVISTIC STABLE WHITE DWARFS 

A. Inverse /3-decay instability 

It is known that white dwarfs may become unstable 
against the inverse /3-decay process [Z, A) ^ [Z — 1, A) 
through the capture of energetic electrons (see e.g. |38l - 
l4l|V In order to trigger such a process, the electron 
Fermi energy must be larger than the mass difference 
between the initial nucleus (Z, A) and the final nucleus 
{Z — \,A). We denote this threshold energy as e§. Usu- 
ally it is satisfied e§_2 < e§ and therefore the initial 
nucleus undergoes two successive decaro, i.e. [Z, A) — > 
[Z -l,A) ^ [Z - 2, A) (see e.g. [1, |^). Some of the 
possible decay channels in white dwarfs with the corre- 
sponding known experimental threshold energies e§ are 
listed in Table [ill The electrons in the white dwarf may 
eventually reach the threshold energy to trigger a given 
decay at some critical density /cfj-jj. Configurations with 
p > become unstable (see [1, Hlj for details) . 

Within the uniform approximation, e.g. in the case of 
the Salpeter equation of state ^ , the critical density for 
the onset of inverse /3-decay is given by 



/iws, we recover the for the density at which the electron Fermi energy ([5; 



/3,unif 



Ar M„, 



-[(4)^-K2m.c^e§]3/^ 



„^,unif gjygjj ^y. ([gg|) and in the relativistic Feynman- 
^/3,roiFMT seen that 



(90) 



equals e^. 

In Tabic [III we show, correspondingly to each threshold 
energy e^, the critical density both in the Salpeter case 

Pent 

Metropolis- Teller case p'^j.l^' 

Pcrlt^^^'^ > pfrit"'^ as one should expect from the fact 
that, for a given density, the electron density at the 
Wigner-Seitz cell boundary satisfies n™'^'^"^ < n™'^. 
This means that, in order to reach a given energy, 
the electrons within the relativistic Feynman- Metropolis- 
Teller approach must be subjected to a larger den- 
sity with respect to the one given by the approximated 
Salpeter analytic formula (pi)). 



B. General relativistic instability 

The concept of the critical mass has played a major 
role in the theory of stellar evolution. For Newtonian 
white dwarfs the critical mass is reached asymptotically 
at infinite central densities of the object. One of the 
most important general relativistic effects is to shift this 
critical point to some finite density p^,^- 

This general relativistic effect is an additional source 
of instability with respect to the already discussed insta- 
bility due to the onset of inverse /3-decay which, contrary 
to the present general relativistic one, applies also in the 
Newtonian case by shifting the maximum mass of New- 
tonian white dwarfs to finite densities (see e.g. {4]|). 



where Eq. ([15)) has been used. 

Because the computation of the electron Fermi en- 
ergy within the relativistic Feynman-Metropolis- Teller 
approach \it\ involves the numerical integration of the 
relativistic Thomas-Fermi equation (|55p. no analytic ex- 
pression for can be found in this case. The critical 
density pfr't'^'^'^"'" is then obtained numerically by looking 



C. Numerical results 

In Figs. [5[[TU[we have plotted the mass-central density 
relation and the mass-radius relation of general relativis- 
tic ^He, ^^C, ^^O and ^^Fe white dwarfs. In particu- 
lar, we show the results for the Newtonian white dwarfs 
of Hamada and Salpeter [ll|, for the Newtonian white 
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H&S 
Pent 




PMTrol 
Pcrit 






1.37 X 10" 


1.44064 


1.56 X 10" 


1.40906 


12^ 


3.88 X 10" 


1.41745 


2.12 X 10" 


1.38603 


16o 


1.89 X 10" 


1.40696 


1.94 X 10" 


1.38024 




1.14 X 10^ 


1.11765 


1.18 X 10^ 


1.10618 



TABLE III: Critical density and corresponding critical mass 
for the onset of gravitational collapse of the Newtonian *He, 
^^C, ^^O and ^®Fe white dwarfs of Hamada [l^, based on the 
Salpeter equation of state :9j, and of the corresponding gen- 
eral relativistic configurations obtained in this work based on 
the relativistic Feynman-Metropolis- Teller equation of state 
[l3|. Densities are in g/cm'^ and masses in solar masses. For 
the sake of comparison, the critical mass of Stoner ([T]) and of 
the one of Chandrasekhar-Landau ^ are M^^,*°"" ~ 1.72Mq 
and M^^!;-^ ~ 1.45A/0, for the average molecular weight 
H = Ar/Z = 2. 



dwarfs of Chandrasekhar [3[ and the general relativistic 
configurations obtained in this work based on the rela- 
tivistic Feynman-Metropolis- Teller equation of state [l3| ■ 
Since our approach takes into account sclf-consistently 
both /3-decay equilibrium and general relativity, we can 
determine if the critical mass is reached due either to 
inverse /3-decay instability or to the general relativistic 
instability. 

A comparison of the numerical value of the critical 
mass as given by Stoner f3|, Eq. ((T)), by Chandrasekhar 
[1| and Landau 6], Eq. by Hamada and Salpeter [11] 
and, by the treatment presented here can be found in 
Table inn 

From the numerical integrations we have obtained: 

1. ^He and ^^C white dwarfs satisfy p^?^ < pf^j^ (see 
Figs. [3Hni and Tables HIl and IIIll) . so they are unsta- 
ble with respect to general relativistic effects. The 
critical density of ^^C white dwarfs is ~ 2.12 x 10^° 
g/cm'^, to be compared with the value 2.65 x 10^" 
g/cm'^ obtained from calculations based on general 
relativistic corrections to the theory of polytropes 
(see e.g. gg). 

2. White dwarfs composed of heavier material than 
"'^^C, e.g. ^^O and ^^Fe are unstable due to inverse 
/9-decay of the nuclei (see Figs. [THTOI and Tables HIl 
and mil) . It is worth to notice that the correct eval- 
uation of general relativistic effects and of the com- 
bined contribution of the electrons to the energy- 
density of the system introduce, for ^^C white 
dwarfs, a critical mass not due to the inverse beta 
decay. When the contribution of the electrons to 
the energy-density is neglected ( e.g. Chandrasekhar 
[1| and Hamada and Salpeter [ll|, see Eq. ((T5|) ) 
the critical density for Carbon white dwarfs is de- 
termined by inverse beta decay irrespective of the 
effects of general relativity. 



3. It can be seen from Figs. ISHlOlthat the drastic de- 
crease of the Salpeter pressure at low densities (see 
d, [13] and Table H] for details) produces an under- 
estimate of the mass and the radius of low density 
(low mass) white dwarfs. 

4. The Coulomb effects are much more pronounced in 
the case of white dwarfs with heavy nuclear com- 
positions e.g. ^^Fe (see Figs. 1^ and ITU)). 



V. CONCLUSIONS 

We have addressed the theoretical physics aspects of 
the white dwarf configurations of equilibrium, quite apart 
from the astrophysical application. 

The recently accomplished description of a compressed 
atom within the global approach of the relativistic Feyn- 
man. Metropolis and Teller [l^l has been here solved 
within the Wigner-Seitz cell and applied to the construc- 
tion of white dwarfs in the framework of general rela- 
tivity. From a theoretical physics point of view, this is 
the first unified approach of white dwarfs taking into ac- 
count consistently the gravitational, the weak, the strong 
and the electromagnetic interactions, and it answers open 
theoretical physics issues in this matter. No analytic for- 
mula for the critical mass of white dwarfs can be derived 
and, on the contrary, the critical mass can obtained only 
through the numerical integration of the general relativis- 
tic equations of equilibrium together with the relativistic 
Feynman-Metropolis- Teller equation of state. 

The value of the critical mass and the radius of white 
dwarfs in our treatment and in the Hamada and Salpeter 
PH treatment becomes a function of the composition of 
the star. Specific examples have been given in the case of 
white dwarfs composed of "^He, ^^C, ^^O and ^^Fe. The 
results of Chandrasekhar, of Hamada and Salpeter and 
ours have been compared and contrasted (see Table IIIII 
and Figs.EHini)- 

The critical mass is a decreasing function of Z and 
Coulomb effects are more important for heavy nuclear 
compositions. The validity of the Salpeter approximate 
formulas increases also with Z, namely for heavy nuclear 
compositions the numerical values of the masses as well as 
of the radii of white dwarfs obtained using the Salpeter 
equation of state are closer to the ones obtained from 
the full numerical integration of the general relativistic 
treatment presented here. 

Turning now to astrophysics, the critical mass of white 
dwarfs is today acquiring a renewed interest in view of its 
central role in the explanation of the supernova phenom- 
ena [isl - lisi . The central role of the critical mass of white 
dwarfs as related to supernova was presented by F. Hoyle 
and W. A. Fowler [4^ explaining the difference between 
type I and type II Supernova. This field has developed in 
the intervening years to a topic of high precision research 
in astrophysics and, very likely, both the relativistic and 
the Coulomb effects outlined in this article will become 
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FIG. 3: Mass in solar masses as a function of the central density in the range (left panel) 10^-10* g/cm'^ and in the range (right 
panel) 10^-5 x lO" g/cm^ for "He white dwarfs. The solid curve corresponds to the present work, the dotted curves are the 
Newtonian configurations of Hamada and Salpeter and the dashed curve are the Newtonian configurations of Chandrasekhar. 




FIG. 4: Mass in solar masses as a function of the radius in units of 10* km for *He white dwarfs. The left and right panels show 
the configurations for the same range of central densities of the corresponding panels of Fig. O The solid curve corresponds to 
the present work, the dotted curves are the Newtonian configurations of Hamada and Salpeter and the dashed curve are the 
Newtonian configurations of Chandrasekhar. 



topic of active confrontation between theory and obser- 
vation. For instance, the underestimate of the mass and 
the radius of low density white dwarfs within the Hamada 
and Salpeter treatment fill (see Figs. [5 HTU|) leads to the 
possibility of a direct confrontation with observations in 
the case of low mass white dwarfs e.g. the companion of 
the Pulsar Jl 141-6545 

We have finally obtained a general formula in Eq. ((76)) 
as a "first integral" of the general relativistic equations 
of equilibrium. This formula relates the chemical po- 
tential of the Wigner-Seitz cells, duly obtained from the 
relativistic Feynman-Metropolis- Teller model [l3| taking 
into account weak, nuclear and electromagnetic interac- 
tions, to the general relativistic gravitational potential 
at each point of the configuration. Besides its esthetic 



value, this is an important tool to examine the radial de- 
pendence of the white dwarf properties and it can be also 
applied to the crust of a neutron star as it approaches to 
the physical important regime of neutron star cores. 

The formalism we have introduced allows in principle 
to evaluate subtle effects of a nuclear density distribu- 
tion as a function of the radius and of the Fermi en- 
ergy of the electrons and of the varying depth of the 
general relativistic gravitational potential. The theoreti- 
cal base presented in this article establishes also the cor- 
rect framework for the formulation of the more general 
case when finite temperatures and magnetic fields are 
present. This treatment naturally opens the way to a 
more precise description of the crust of neutron stars, 
which will certainly become an active topic of research 



14 




FIG. 5: Mass in solar masses as a function of the central density in the range (left panel) 10^-10* g/cm'^ and in the range 
(right panel) 10*-10" g/cm^ for ^^C white dwarfs. The solid curve corresponds to the present work, the dotted curves are the 
Newtonian configurations of Hamada and Salpeter and the dashed curve are the Newtonian configurations of Chandrasekhar. 




FIG. 6: Mass in solar masses as a function of the radius in units of 10* km for ^^C white dwarfs. The left and right panels show 
the configurations for the same range of central densities of the corresponding panels of Fig. [5] The solid curve corresponds to 
the present work, the dotted curves are the Newtonian configurations of Hamada and Salpeter and the dashed curve are the 
Newtonian configurations of Chandrasekhar. 



in view of the recent results by S. Goriely et al. [5l|, |52| 
and by J. M. Pearson et al. [5J| on the importance of the 
Coulomb effects in the r-process nucleosynthesis of the 
crust material during its post-ejection evolution in the 
process of gravitational collapse and/or in the merging 
of neutron star binaries. 



Acknowledgments 



We would like to thank Professor D. Arnctt for discus- 
sions. 



[1] R. H. Fowler, Mon. Not. Roy. Astr. Soc. 87, 114 (1926). 
[2] E. C. Stoner, Philosophical Magazine (Series 7) 7, 63 
(1929). 

[3] S. Chandrasekhar, Astrophys. J. 74, 81 (1931). 

[4] E. A. Milne, Mon. Not. Roy. Astr. Soc. 91, 4 (1930). 



[5] R. Emden, Gaskugeln Anwendungen der Mechamschen 
Wdrmetheorie auf Kosmologische und Meteorologische 
Probleme (Leipzig, Teubner, Berlin, 1907). 

[6] L. D. Landau, Phys. Z. Sowjetunion 1, 285 (1932). 

[7] K. C. Wah, Physics Today 35, 33 (1982). 



15 




FIG. 7: Mass in solar masses as a function of the central density in the range (left panel) 10^-10* g/cm^ and in the range 
(right panel) 10*-10" g/cm^ for ^''O white dwarfs. The solid curve corresponds to the present work, the dotted curves are the 
Newtonian configurations of Hamada and Salpeter and the dashed curve are the Newtonian configurations of Chandrasekhar. 




FIG. 8: Mass in solar masses as a function of the radius in units of 10* km for ^^O white dwarfs. The left and right panels show 
the configurations for the same range of central densities of the corresponding panels of Fig. [T] The solid curve corresponds to 
the present work, the dotted curves are the Newtonian configurations of Hamada and Salpeter and the dashed curve are the 
Newtonian configurations of Chandrasekhar. 



[8] D. Boccaletti and R. RufEni, Fermi and Astrophysics 

(World Scientific, Singapore, 2010). 
[9] E. E. Salpeter, Astrophys. J. 134, 669 (1961). 
[10] Y. I. Frenkel, Zeit. fur Phys. 50, 234 (1928). 
[11] T. Hamada and E. E. Salpeter, Astrophys. J. 134, 683 
(1961). 

[12] E. Olson and M. Bailyn, Phys. Rev. D 12, 3030 (1975). 
[13] E. Olson and M. Bailyn, Phys. Rev. D 13, 2204 (1976). 
[14] M. Rotondo, J. A. Rueda, R. Ruffini, and S. Xue, Phys. 

Lett. B 701, 667 (2011). 
[15] R. Ruffini, in Exploring the universe: a Festschnft 

m honor of Riccardo Giaccom, edited by H. Gursky, 

R. Ruffini, & L. Stella (2000). 
[16] G. Bertone and R. Ruffini, Nuovo Cimento B Serie 115, 

935 (2000). 

[17] M. Rotondo, J. A. Rueda, R. Ruflini, and S.-S. Xue, 



Phys. Rev. C 83, 045805 (2011), 0911.4622. 
[18] R. P. Feynman, N. Metropolis, and E. Teller, Phys. Rev. 

75, 1561 (1949). 
[19] G. Baym, C. Pethick, and P. Sutherland, Astrophys. J. 

170, 299 (1971). 
[20] L. D. Landau and E. M. Lifshitz, Statistical phystcs. 

Parti (Pergamon Press, Oxford, 1980). 
[21] J. Ferreirinho, R. Ruffini, and L. Stella, Phys. Lett. B 

91, 314 (1980). 
[22] R. Ruffini and L. Stella, Phys. Lett. B 102, 442 (1981). 
[23] P. A. M. Dirac, Proc. Cambridge Phil. Soc. 26, 376 

(1930). 

[24] A. B. Migdal, V. S. Popov, and D. N. Voskresenskii, So- 
viet Journal of Experimental and Theoretical Physics 45, 
436 (1977). 

[25] J. C. Slater and H. M. Krutter, Phys. Rev. 47, 559 



16 




FIG. 9: Mass in solar masses as a function of the central density in the range (left panel) 10^-10* g/cm'^ and in the range (right 
panel) 10*-3 x lO'' g/cm^ for ^"^Fe white dwarfs. The sohd curve corresponds to the present work, the dotted curves are the 
Newtonian configurations of Hamada and Salpeter and the dashed curve are the Newtonian configurations of Chandrasekhar. 




FIG. 10: Mass in solar masses as a function of the radius in units of 10* km for ^^Fe white dwarfs. The left and right panels show 
the configurations for the same range of central densities of the corresponding panels of Fig. [S] The solid curve corresponds to 
the present work, the dotted curves are the Newtonian configurations of Hamada and Salpeter and the dashed curve are the 
Newtonian configurations of Chandrasekhar. 



(1935). 

[26] R. Ruffini, in Path Integrals - New Trends and Perspec- 
tives (2008), pp. 207-218. 

[27] V. S. Popov, M. Rotondo, R. Ruflini, and S. Xue, in Pro- 
ceedings of the 1st Galileo-Xu Guantqi Meettng (2009), 
Int. J. Mod. Phys. Conf. S. in press, edited by D. Blair, 
R. Ruffini, J. Yipeng and S.-S. Xue (2011), arXiviastro- 
ph/0903.3727. 

[28] R. C. Tolman, Physical Review 55, 364 (1939). 

[29] J. R. Oppenheimer and G. M. VolkoflF, Phys. Rev. 55, 
374 (1939). 

[30] O. Klein, Reviews of Modern Physics 21, 531 (1949). 
[31] R. C. Tolman, Physical Review 35, 904 (1930). 
[32] R. C. Tolman and P. Ehrenfest, Physical Review 36, 1791 
(1930). 

[33] S. Chandrasekhar, Mon. Not. Roy. Astr. Soc. 91, 456 



(1931). 

[34] S. Chandrasekhar, Mon. Not. Roy. Astr. Soc. 95, 207 
(1935). 

[35] S. Chandrasekhar, An introduction to the study of stellar 

structure (1939). 
[36] A. S. Eddington, Sir, Mon. Not. Roy. Astr. Soc. 95, 194 

(1935). 

[37] I. Ciufolini and R. Ruffini, Astrophys. J. 275, 867 (1983). 
[38] F. Hund, Erg. d. exacten Natwis. 15, 189 (1936). 
[39] L. D. Landau, Nature 19, 333 (1938). 
[40] I. B. Zel'Dovich, Soviet Journal of Experimental and 

Theoretical Physics 6, 760 (1958). 
[41] B. K. Harrison, M. Wakano, and J. A. Wheeler, Onzieme 

Conseil de Physisque de Solvay (1958). 
[42] S. L. Shapiro and S. A. Teukolsky, Black holes, white 

dwarfs, and neutron stars: The physics of compact oh- 



17 



jects (1983). 

[43] G. Audi, A. H. Wapstra, and C. Thibault, Nuclear 
Physics A 729, 337 (2003). 

[44] A. H. Wapstra and K. Bos, Atomic Data and Nuclear 
Data Tables 19, 175 (1977). 

[45] M. M. Phillips, Astrophys. J. Lett. 413, L105 (1993). 

[46] A. G. Riess, A. V. Filippenko, P. Challis, A. Clocchiatti, 
A. Diercks, P. M. Garnavich, R. L. Gilliland, C. J. Hogan, 
S. Jha, R. P. Kirshner, et al.. Astronomical Journal 116, 
1009 (1998), arXiv:astro-ph/9805201. 

[47] S. Perlmutter, G. Aldering, G. Goldhaber, R. A. Knop, 
P. Nugent, P. G. Castro, S. Deustua, S. Fabbro, A. Goo- 
bar, D. E. Groom, et al., Astrophys. J. 517, 565 (1999), 
arXiv:astro-ph/9812133. 

[48] A. G. Riess, L. Strolger, J. Tonry, S. Casertano, H. C. 
Ferguson, B. Mobasher, P. Challis, A. V. Filippenko, 
S. Jha, W. Li, et al., Astrophys. J. 607, 665 (2004), 
arXiv:astro-ph/0402512. 

[49] F. Hoyle and W. A. Fowler, Astrophys. J. 132, 565 
(1960). 

[50] M. Kramer, private communication (2010). 

[51] S. Goriely, N. Chamel, H.-T. Janka, and J. M. Pearson, 

A. & A. 531, A78+ (2011), 1105.2453. 
[52] S. Goriely, A. Bauswein, and H. -Thomas Janka, ArXiv 

e-prints (2011), 1107.0899. 



[53] J. M. Pearson, S. Goriely, and N. Chamel, Phys. Rev. C 

83, 065810 (2011). 
[54] M. Nauenberg, Journal for the History of Astronomy 39, 

297 (2008). 

[55] J. L. Heilbron, Historical Studies in the Physical Sciences 
13, 261 (1983). 

[56] In doing this, Stoner used what later became known as 
the exclusion principle, generally attributed in literature 
to Wolfgang Pauli. For a lucid and scientifically correct 
historical reconstruction of the contributions to the criti- 
cal mass concept see [Hj . For historical details about the 
exclusion principle see also [sst . 

[57] At the time a 20 years old graduate student coming to 
Cambridge from India. 

[58] The dispute reached such a heated level that Chan- 
drasekhar was confronted with the option either to 
change field of research or to leave Cambridge. As is well 
known he chose the second option transferring to Yerkes 
Observatory near Chicago where he published his results 
in his classic book [ssl ]. 

[59] It goes to Eddington credit, at the time Plumian Pro- 
fessor at Cambridge, to have allowed the publication of 
the Chandrasekhar work although preceded by his own 
critical considerations ,36] . 



